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The existence of multidimensional lattice compactons in the discrete nonlinear Schrodinger equa¬ 
tion in the presence of fast periodic time modulations of the nonlinearity is demonstrated. By 
averaging over the period of the fast modulations, a new effective averaged dynamical equation 
arises with coupling constants involving Bessel functions of the first and zeroth kind. These terms 
allow one to solve, at this averaged level, for exact discrete compacton solution configurations in 
the corresponding stationary equation. We focus on seven types of compacton solutions: single site 
and vortex solutions are found to be always stable in the parametric regimes we examined. Other 
solutions such as double site in- and out-of-phase, four site symmetric and anti-symmetric, and a 
five site compacton solution are found to have regions of stability and instability in two-dimensional 
parametric planes, involving variations of the strength of the coupling and of the nonlinearity. We 
also explore the time evolution of the solutions and compare the dynamics according to the averaged 
with those of the original dynamical equations without the averaging. Possible observation of com¬ 
pactons in the BEG loaded in a deep two-dimensional optical lattice with interactions modulated 
periodically in time is discussed. 

PACS numbers: 42.65.-k, 42.81.Dp, 03.75.Lm 


I. INTRODUCTION 


Periodic management of parameters of nonlinear lat¬ 
tices (and continua) is a very attractive technique for 
the generation of new types of systems or excitations 
with interesting localization properties [T]. For exam¬ 
ple, dynamic suppression of (atom and light, respec¬ 
tively) interwell tunneling has been reported experimen¬ 
tally in Bose-Einstein condensates [2] and in optical 
waveguide arrays [3]. On the other hand, very fast pe¬ 
riodic time variation of the nonlinearity, also called the 
strong nonlinearity management (SNLM) technique, has 
been recently shown to be quite effective towards induc¬ 
ing compactly supported (so-called compacton) solutions 
in lower dimensional discrete systems such as one dimen¬ 
sional (ID) one and two component Bose-Einstein con¬ 
densates (BECs) in optical lattices or arrays of nonlinear 
optical waveguides described by the discrete nonlinear 
Schrodinger (DNLS) equation [H [5]. The main feature 
of such solutions, unlike other types of nonlinear exci¬ 
tations such as discrete breathers and intrinsic localized 
modes, is that their amplitudes decay to zero sharply 
without any tails. It has been shown for the ID case 
that the lack of exponential tails is a consequence of the 
nonlinear dispersive interaction induced by the SNLM 
which permits the vanishing of the intersite tunneling at 
compacton edges. The absence of tails imply that com¬ 
pactons cannot interact with each other until they are in 


contact. In the lattice realm, this potentially allows for 
the possibility of “maximal localization” in the form of a 
single site solution. 

Compactons are quite generic in nonlinear continua [6l 
|7] (in one and also higher [8] dimensions) and lat¬ 
tices [9l [To] (including compact breathers m) bearing 
nonlinear dispersion. The difficulty of implementing this 
condition in physical contexts has largely restricted the 
investigations mainly to the mathematical side, although 
this situation is now rapidly evolving. Besides Bose- 
Einstein condensates under SNLM, physical compactons 
were suggested to appear also in exciton-polariton con¬ 
densates [T2| and nearly-compact (doubly exponential) 
traveling waves have also been identified in the realm 
of granular crystals bearing purely nonlinear interac¬ 
tions m- Finally, another area of significant interest 
has recently arisen where special solutions in the form 
of discrete compactons may emerge. This is due to the 
existence of so-called fiat bands in the linear dispersion 
relation (due to the geometric characteristics of the cor¬ 
responding lattice, such as e.g. the Kagome lattice) [14] . 
A realization of this type emerged very recently in the 
realm of the so-called Lieb photonic lattices m- 

Among the above different (Klein-Gordon, nonlinear 
Schrodinger, Fermi-Pasta-Ulam, etc.) model equations 
in which compactly supported structures have been pro¬ 
posed DNLS ones are, arguably, of the most widely ap¬ 
plicable as models both for BEGs and for nonlinear op¬ 
tics. In the presence of SNLM, this type of system sup- 
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ports discrete compactons via a site-dependent nonlinear 
rescaling of the inter-well tunneling (see [4] for details). A 
similar approach, applied to the quantum version of the 
DNLS model, i.e., to the Bose-Hubbard model with time 
dependent onsite interaction, has been shown to be quite 
effective for creating new quantum phases in BECs with 
compacton-type excitations [ISIIIT] as well as for the gen¬ 
eration of density dependent synthetic gauge fields m- 

These studies have mainly focused on the one di¬ 
mensional case. On the other hand, it is known that 
compactons can exist also in multidimensional con¬ 
texts. In particular, multidimensional compactons have 
been investigated in continuous models such as (two- 
dimensional) variants of the Korteweg-de Vries equations 
with nonlinear dispersion [ 8 ], relativistic scalar field the¬ 
ories in two dimensions m, Klein-Gordon and nonlin¬ 
ear Schrodinger equations with sub-linear forces, etc.. In 
these last cases, however, the mechanism leading to com- 
pacton formation is not the nonlinear dispersion but the 
presence of a sub-quadratic interaction potential which 
enforces compact patterns with sharp fronts [201 El]. In 
this context, existence of vortex compactons was also 
demonstrated to be possible, although only with a finite 
lifetime [ 21 ]. 

In the discrete case, multidimensional compactons 
have been scarcely investigated and only in the special 
case example of compact coherent structures in the pres¬ 
ence of a flat band of the linear spectrum nms]. To 
the best of our knowledge, general case examples (dif¬ 
ferent numbers of sites, including ones bearing vorticity) 
of multidimensional discrete compactons of the nonlin¬ 
ear Schrodinger type, amenable to physical applications 
in BECs and nonlinear optics, have not yet been investi¬ 
gated. 

The aim of this paper is to demonstrate and system¬ 
atically explore the existence and stability of multidi¬ 
mensional lattice compactons in the discrete nonlinear 
Schrodinger equation in the presence of periodic time 
modulations of the nonlinearity i.e., in a 2 -dimensional 
generalization of |4] bearing the potential for a wide range 
of additional structures. In this work we concentrate 
mainly on compactons localized on no more than five 
interacting neighboring sites of a two-dimensional (2D) 
square lattice; of particular interest, in addition to the 
stable single site compacton, is the stable four site com- 
pacton vortex solution i.e., bearing a vortical phase struc¬ 
ture. In particular, the existence and stability properties 
of DNLS compacton excitations are investigated both 
for generic time dependent nonlinear modulations and 
in the SNLM limit for which an effective averaged DNLS 
model is derived. We show that single site compactons 
are stable in the whole parameter space, however, two site 
compactons have finite stability ranges (although bearing 
some nontrivial differences from the “standard” DNLS 
model), different for symmetric and anti-symmetric types 
and with different dynamical features. Interestingly, sta¬ 
tionary three site compactons with real amplitude cannot 
exist in the square geometry, while four site compactons 


of the symmetric or anti-symmetric types exist but have 
finite regions of instability similar to the double site con¬ 
figurations. On four nearest-neighboring sites, however, 
for all values in the parameter space one has a stable 
vortex compacton in which the phase of the wavefunc- 
tion increases by 7r/2 moving clockwise from one corner 
to the next of the square. Since compactons interact only 
when they are in contact, one can obviously construct ar¬ 
bitrary single and two site compacton patterns by placing 
them on non interacting sites (such as for example next- 
neighboring sites along a diagonal). We also show that 
five site compactons with C 4 symmetry can exist and can 
be stable; their regions of instability are finite, similar to 
other configurations’ instability. 

The paper is organized as follows. In section II we in¬ 
troduce the model equation of a 2 D DNLS with SNLM 
and discuss the theoretical derivation of the effective av¬ 
eraged equations with nonlinear dispersion. In section III 
we derive the existence conditions for exact compacton 
solutions of the averaged DNLS equation and in section 
IV we study numerically the linear stability properties 
and compare results with direct numerical integrations of 
the original (unaveraged) DNLS system. In the last sec¬ 
tion we discuss briefly the potential future experimental 
implementations, summarize our main results and pro¬ 
vide some pointers for future work. 


II. MODEL AND THEORETICAL SETUP 

Consider the following 2 D DNLS equation [ 22 ] 

T ^("^n+ljm T "^71 — 1,771) T ^("^ 77 , 771+1 T "^77,771—1) 

T('70 T 1^77,777 I '^77,777 — O5 (l) 

which serves as a model for the dynamics of BEG in opti¬ 
cal lattices subjected to SNLM (through varying the in¬ 
teratomic scattering length by external time dependent 
magnetic fields via a Eeshbach resonance) [23] , as well as 
for light propagation in 2D optical waveguide arrays [24] . 
In the latter case, where this type of modulation has been 
realized not only in discrete but also in continuum me¬ 
dia (in both one and higher dimensions) [25], the evo¬ 
lution variable is the propagation distance. Hence, here 
the SNLM consists of periodic space variations of the 
Kerr nonlinearity along the propagation direction with 
the coupling constants k, k quantifying the tunneling be¬ 
tween adjacent sites along n and m directions, respec¬ 
tively, 7 o denoting the onsite constant nonlinearity and 
7 (t) representing the time dependent modulation (of the 
interatomic interactions or of the refractive index, in 
atomic and optical settings, respectively). In the fol¬ 
lowing we assume a strong management case with 7 (t) 
being a periodic, e.g. + Tq), and rapidly 

varying function. As a prototypical example, we use 
7 (t) = ^cos(Ur), with 71 ^ 0 ( 1 ), 5 <C 1 , r = t/e 
denoting the fast time variable and T = 27 t/Q the pe¬ 
riod with respect to r (Tq = sT). In the following we 
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take, for simplicity, the coupling constants k^k to he the 
same for the two directions (assuming square symmetry): 
k = k = n. 

The existence of compacton solutions of Eq. 0 in 
the SNLM limit can be inferred from (and analyzed in 
the context of) an effective averaged 2 D DNLS equation 
obtained by averaging out the fast time r. To that ef¬ 
fect, it is convenient to perform the transformation [26] 
Un,m{t) = with T = dt j{t) = 

sin (Dr), which allows one to rewrite Eq. Q as 

~ I )t 'To|'^n,m| "^71,1711 ( 2 ) 

with 

, and (3) 

On the other hand, {\Vn,m\‘^)t = {Vn,mVn,m + '>^n,mVn,m) = 
with the Star denoting complex 
conjugation. Substituting this expression into Eq. 0 
and averaging the resulting equation over the period T 
of the rapid modulation, we obtain 

^hn,m = m|nn,mp(rx) - -/^(X) 

'To|'^n,m| (4) 

with (•) = ^ (')dr denoting the fast time average. 

The averaged terms in Eq. 0 can be calculated by 
means of the elementary integrals = Jo(o^^±), 

_ ±iaJi(ad±), with Ji being Bessel functions 
of order i = 0,1 and a = 71 /D, thus giving 

-^('^) (^) 

with 

F{v) = - 7 o|'^n,mP'^n,m 

[('^n+l,m'^n,m ~h '^n-\-l,rn^'ri,rn) ) 

+ ('^n,m+l'^n,m + '^n )Ji{a0i) 

+ ('^n —~h '^n—l,m^n,m)^l (<^^— ) 

~h{Vn,m-l'^ri,m ~h )] (^) 

^['^n+l,m'A)( q^^_I_ ) T '^n,m+l^0(^^n ) 

-hVn-l,mJo{o^O^) + Vn^m-lJo{o^O^ )]. 

Note that the parameters 71 , D ^ 0(1), and the averaged 
equation is valid for times t < 1 /e. This modified DNLS 
equation can be written as ivn,m = ^Hav/^v^ rn^ with the 
averaged Hamiltonian 

Hav = ~ Jo(a^^ ) [nn+l,m'^n,m ~h 

n,m 

-\-I^Jo{(yO^) [Vn^rn-\-l'^n,m ~h '^n,m+l'^ri,m\ + ~^\'^n,m\ }• 

It is interesting to note that this Hamiltonian, except for 
the rescaling of the coupling constants k KJo{aO'^), 
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FIG. 1: (color online) We plot typical examples of exact 2D 
compacton solutions of Eq. 0- In the top left panel we show 
a single site solution, in the bottom left we show a double-site 
in-phase solution, and in the bottom right, we show a double¬ 
site out-of-phase solution. In these three plots colors show 
the (real) values of An,m where blue represents An,m = Aq 
and red An,m = —^o- No color represents zero amplitude. In 
the top right plot we show a vortex solution where the four 
central colored sites all have amplitude Aq, and arg{vn,m) ^ 
{—7r/2, 0,7r/2 , tt} is plotted by the four respective colors: dark 
red (at (n,m) = (0,1)), light red (at (n,m) = (0,0)), light 
blue (at (n,m) = (1,0)), dark blue (at (n,m) = (1,1))- Here 
the length and width of the 2D grid are each L — 7. 


k /^Jo(a 6 >+), is the same as the Hamiltonian of the 
DNLS equation in absence of SNLM (e.g. with 71 = 0). 
A similar rescaling was reported also for the ID case [4]. 

It is also worth noting that while the appearance of the 
Bessel function is intimately connected with harmonic 
modulations, the existence of compacton solutions and 
the lattice tunneling suppression is generic for periodic 
SNLM. 

The generalization of these equations to the case of 3D 
DNLS with cubic lattice symmetry is also quite straight¬ 
forward to derive (omitted here for brevity). In our nu¬ 
merical results below, we will restrict our considerations 
to the numerically more tractable 2D case (also more 
physically realistic at least in the optics realm where the 
^-direction plays the role of the propagation direction). 


III. EXACT 2D COMPACTONS 

To demonstrate the existence of exact compactons in 
the averaged system, we seek stationary solutions of the 
form Vn^m = with An^m ^ C for which Eq. 
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FIG. 2 : (color online) The plots show three additional exam¬ 
ples of exact 2D compacton solutions of Eq. 0 on four and 
five lattice sites. In the left and middle panels symmetric and 
anti-symmetric four site compactons are shown, respectively, 
while the right panel depicts a hve site compacton. Dark blue 
and red colors represent An,m = Aq and An^m = — similar 
to the single and double site plots in Figure The light blue 
color at the center of the right-most plot represents the value 
y for the center site of the hve site solution obtained from Eq. 
(16). Axis labels are the same here as those seen in Figure 


becomes 


Mn,m = F{A) (7) 

for F given by Eq. In the following, we theoretically 
predict and numerically verify that exact solutions of this 
equation can exist in the form of genuine compactons i.e., 
possessing vanishing tails. 

To search for compacton solutions, we begin by apply¬ 
ing 0 at any zero-amplitude site that has exactly one 
out of its four neighbors nonzero. We label this nonzero 
site An,m 7^ 0 and we call it an “edge site” of the com¬ 
pacton. For /i: 7 ^ 0, this gives the condition 

Jo{a\AN,Mf) = 0 ^ \AN,Mf = Zj/a ( 8 ) 

where Zj is the zero of the Bessel function Jq. In other 
words, each edge site of the compacton has amplitude 

Ao ''=■ ^Jzj/a. (9) 

Applying 0 at edge site(s) of the compacton gives addi¬ 
tional condition(s) that depend on the shape of the com¬ 
pacton. We will focus on the following shapes although 
other configurations are possible: a single site, a double 
site in-phase, a double site out-of-phase, three types of 
four site solutions in the form of symmetric, asymmetric 
and vortex excitations, and finally a five site configura¬ 
tion. Figures and show one example of each of these 
configurations. L represents the horizontal (and vertical) 
length of the lattice. 

For a single site compacton solution located at 
(no,mo), this nonzero site is an edge site so it has am¬ 
plitude Aq. Taking ^no,mo = and applying ([^ at 
(no, mo) gives the value 


( 10 ) 


FIG. 3: (color online) Graphical solution of Eq. (16) written 
as G{y) = 0 (left panel) and section at m = 0 of a hve- 
site compacton (right panel) of Eq. 0 for parameter values 
n = 0 . 5,0 = l, 7 o = 2 with amplitudes at edges sites fixed 
in corresponden ce of t he fifth zero, Z 5 , of the Bessel function 
Jo, e.g. Ao = yj z^jOL ~ 3.864. The amplitude of the central 
site is taken as the zero of G(y) at y ~ 3.8695 (black dot in 
the left panel) while the corresponding chemical potential is 
obtained from Eq. (15) as ~ —30.678. 


It is interesting to notice that this solution coincides with 
the one obtained for the one dimensional DNLS under 
SNLM in [4] . This is due to the fact that for a single site 
the compact nature of the solution does not allow it to 
distinguish ID from 2D or even 3D (this is true for a 3D 
cubic lattice as well). Similar to the ID case the single 
site compactons are stable; we will see this in more detail 
in Section m 

For a double-site compacton located at (no,mo) and 
(no + l,mo) (or any other adjacent pair of indices) 
each with amplitude Ao there are the two possibilties: 
Ano+i,mo = i^no,mo where plus denotes an in-phase so¬ 
lution and minus an out-of-phase solution. Using that 
Ji(0) = 0 and Jo(0) = 1, application of 0 at either of 
the two nonzero sites then gives the value 

l_i = zfK- 7oAg (11) 

where the minus corresponds to the in-phase solution and 
the plus to the out-of-phase solution. 

From the above cases one could expect that, although 
with different stability properties, ID compactons lo¬ 
cated on the lattice along n or m directions could be 
also solutions of multidimensional lattices. This how¬ 
ever, except for the cases considered above, is not true 
in general. In this respect, one and two-sites compactons 
are very special, since the number of edge sites in these 
solutions does not change with dimensionality their an¬ 
alytical expressions remain the same in ID, 2D and 3D. 
For compactons with more than two sites, however, this 
is not true. Thus, for three-site compacton with real am¬ 
plitudes (..., 0, Ai, A 2 , Ai, 0,...), in the n direction and 
zero amplitudes in all other sites, for example, one has 
three edge sites in 2D but only two in ID. In contrast 
with the ID case, where solutions with real amplitudes 
exist for arbitrary choices of parameters and are stable 
m, the constraints resulting from Eq. 0 can be satisfied 
only for particular choices of parameters, e.g. real ampli¬ 
tude three-site compactons of 2D DNLS under SNLM are 


II = -7oAo. 
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with the signs plus and minus referring to the symmetric 
and anti-symmetric case, respectively. 

Quite remarkably, the DNLS system under SNLM can 
also support discrete vortex-compactons, an unprece¬ 
dented feature of earlier two-dimensional (continuous or 
discrete) compacton considerations. In particular, for a 
2x2 shaped vortex solution with the top left nonzero 
site located at index (no, uiq), all four sites are edge sites 
with amplitude Aq. Then we can write 

^no,mo ~ ^0? ^no,mQ-\-l — AqC ^ 

^no + l,mo = + + i (13) 

and 0 applied on any of the four vortex sites gives 

H = --foAl. (14) 


FIG. 4: (color online) We plot max(real(z/)) in order to show 
the strength of the instabilities and regions of stability. Pa¬ 
rameter values common across all four plots are e = 0.01 and 
Q = 7 i = 10 so that a = 1. The left column of plots is 
associated with a double site in-phase solution with lattice 
length both horizontally and vertically L = 17. The right 
column of plots is associated with a double site out-of-phase 
solution with lattice length L = 15. The top row corresponds 
to j = 1 and the bottom to j = 2. Notice that the top right 
plot shows more striations in the coloring than the top left 
plot; these further smooth out with higher lattice size L. The 
green dots in the top two plots are parameter values for which 
the unstable solutions are propagated in time; see Figures 
[S] below. 


not generic solutions in parameter space. The situation 
is different for complex and an example of this type will, 
in fact, be given below. 

Genuine 2D compacton solutions first occur on four 
neighboring sites of a lattice cell. These square com- 
pactons can be either symmetric or anti-symmetric (see 
the panels of Fig. [^. In either case all four nonzero sites 
are edge sites with amplitude Aq. Denoting with (no, mo) 
the coordinates of the upper left corner of the four site 
compacton, the symmetric solution has plus signs on the 
four nonzero sites so that 

^no,mo ~ ^no + l,mo ~ ^no,mo + l ~ ^no + l,mo + l “ ^0- 

The anti-symmetric solution has an alternating pattern 
of plusses and minuses of the form 

^no,mo ~ ^no + l,mo + l ~ ^0? 

^no + l,mo ~ ^no,mo + l ~ ^0: 

with vanishing amplitudes, = 0, on all other sites. 

Substituting the above expressions into Eq. 0 one read¬ 
ily gets the following equation for /i to be satisfied 

zb 2 k + ^oTo + /i = 0 ( 12 ) 


With some additional effort it is possible to obtain ex¬ 
act compacton solutions involving more (> 4) nonzero 
sites. As an example, consider the case of a five site 
compacton shaped as in the right panel of Fig. Let 
(no, mo) denote the index of the center site of the con¬ 
figuration and define y = Ano,mo- The four edge sites 
of the five site solution have amplitudes Aq so we take 

^no + l,mo ~ ^no —l,mo ~ ^no,mo + l ~ ^no,mo — l ~ ^0- 

Notice that the four edge sites (dark blue in the right 
panel of Fig. isolate the central (light blue) site from 
the bulk, reproducing the same situation occurring for 
a three-site compacton in ID (in 3D a similar solution 
would imply a seven-site compacton). Plugging the above 
edge site expressions into Eq. 0, applied at any one of 
the four edge sites, then gives an expression fixing the 
chemical potential 


li{y) = -K-^Jo{^{y))-'yoAl + 2aKAoyJi{^{y)), (15) 


where ^{y) = a{A^y‘^). The final constraint is obtained 
by applying Eq. Q at the center site (no, mo). Combin¬ 
ing the result with the expression in Eq. (15) for /i shows 
that the amplitude y at the center site must be a solution 
of the equation G(^) = 0 for 


G{y) ^ Jo(e(y)) - 

\ y AoJ a 

+ 10 aKAoyJi{^{y)). 


(16) 


Although Eq. (16) does not yield a simple analytical ex¬ 


pression for it can be easily solved numerically provid¬ 
ing an exact five-site compacton with chemical potential 
given by Eq. (15). In general, for fixed parameters and 


a fixed value for Aq = more than one zero can 

exist for G (see left panel of Eig. [^, although not all of 
them necessarily lead to stable compactons. It is remark¬ 
able, however, that some of them can have non-vanishing 
intervals/ranges of stability in parameter space. We now 
turn to the details of the stability considerations. 
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FIG. 5: (color online) The panels show eigenvalues in the 
complex plane for selected parameters corresponding to a hor¬ 
izontal cut on the top left panel of Figure (in-phase solu¬ 
tion), obtained by setting 70 = 1.75. For ac = 0 (not shown) 
the eigenvalues are purely imaginary and four eigenvalues are 
zero. For small nonzero k there are only two zero eigenval¬ 
ues and the other two previously zero eigenvalues have begun 
to increase in magnitude on the real axis; this is shown for 
K, = 0.3 in the top left panel. As k increases the solution 
stabilizes and these real eigenvalues move back towards zero 
horizontally along the real axis. For even higher k values (at 
which the solution remains stable), the two eigenvalues in¬ 
crease in magnitude again from zero, but this time along the 
imaginary axis; this is shown for k, = 0.5 in the upper right 
panel. These eigenvalues merge with the band of eigenvalues 
on the imaginary axis meanwhile the (continuous spectrum) 
band lengthens; this is shown for = 1 in the bottom left 
panel. At even higher k values, two complex conjugate purely 
imaginary eigenvalues emerge from the band at a magnitude 
higher than any eigenvalue in the band; this is depicted in the 
lower right panel for = 5. The eigenvalues in the bottom 
right panel that appear to have very small nonzero real part 
are seen to actually have zero real part when the eigenvalues 
are computed for a much higher lattice length L. 


FIG. 6 : (color online) The panels show the eigenvalues in the 
complex plane for selected parameters corresponding to a hor¬ 
izontal cut on the top right panel of Figure (out-of-phase 
solution), obtained by setting 70 = 3. For At = 0 (not shown) 
the eigenvalues are purely imaginary and four eigenvalues are 
zero. For small nonzero At there are only two zero eigenvalues 
and the other two previously zero eigenvalues have begun to 
increase in magnitude on the imaginary axis; this is shown 
for At = 0.35 in the top left panel. After the two merge with 
the (continuous spectrum) band of eigenvalues on the imagi¬ 
nary axis a quartet of complex eigenvalues then emerges for 
higher At; this is depicted for At = 0.45 in the top right panel. 
For higher At values the four complex eigenvalues move in the 
complex plane, first increasing in real and imaginary parts 
then decreasing in real part until they merge again with the 
imaginary axis (meanwhile the continuous spectrum band on 
the imaginary axis lengthens); the bottom two plots show this 
restabilization with At = 1.38 in the bottom left and At = 1.45 
in the bottom right. At even higher At values (at which the 
solution remains stable) two complex conjugate purely imagi¬ 
nary eigenvalues emerge from the band at a magnitude higher 
than any eigenvalue in the band, similar to the lower right 
panel of Figure 


IV. STABILITY ANALYSIS AND NUMERICAL 
SIMULATIONS 


In order to examine stability of the solutions found in 
Section EH we set 


= {An,m + 6 


(17) 


where both An^rn and ^n,m are independent of time t. 
One can then show that for '0n,m = cin,m the real 


vector 


(with 2 L^ components) is an eigenvector 


of the 2L? X 2LP‘ matrix 


M = 


da 

dRejFj'iP)) 

da 


dIm{F(2p)) 

db 

dRe{Fi'iP)) 

db 


(18) 


with eigenvalue u. 

We find that the single site and vortex solutions 
are stable for all parameter values Xq so that 
max(real(z/)) = 0; in other words, the nonzero sites’ time 
evolution plots according to the averaged Eq. ^ even 
under perturbation does not lead to growth. The time 
evolution of the single site or vortex solution according to 
the original Eq. 0 with the time-periodic nonlinearity 
gives the expected result of oscillations about the stable 
solution given by the averaged equation. The double site 
solutions have regions of instability which depend on the 
system’s parameters. Eigurej^shows the maximal growth 
rate of the instability as a function of the dc strength of 
the nonlinearity and of the coupling constant. We see 
that the in-phase solutions present instability intervals 
near the uncoupled limit (but contrary to what is the 
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FIG. 7: (color online) The plots show the result of propaga¬ 
tion of the double site in-phase exact compacton solution in 
the time parameter t according to each of the equations Q 
and 0 . The values of e, Q, 71 , a are the same as in Figure^ 
these values give Aq ^ 1.5507 for j — 1 through Eq. 0 . In 
the bottom left panel the plot shows log 3 ^Q(|?Xn,m( 0 )|) plotted 
as a function of n^m where r^n,m( 0 ) = 'Cn,m( 0 ) is the initial 
double site in-phase compacton configuration at t = 0 ; this 
plot, of course, shows an order 0 logarithmic amplitude at the 
two excited sites and vanishing amplitude elsewhere. The in¬ 
stability strength for this solution is max(real(z/)) ~ 0.5327. 
In the bottom right panel the plot shows log^Q(|r'n,m(120)|) 
where Vn,m{t) is determined from r’n,m( 0 ) according to the 
averaged equation ([^; the plot of log]^Q(|i/n,m( 120 )|) deter¬ 
mined from the non-autonomous equation 0 is visually in¬ 
distinguishable from the averaged version shown in the bot¬ 
tom right. In the bottom right plot we see small ampli¬ 
tude excitations appearing near the compacton site, and these 
small amplitudes continue to spread out spatially as t in¬ 
creases, while the central sites have amplitudes which ap¬ 
proach nonzero values. The top left panel shows this evolu¬ 
tion of the amplitude of the two excited sites over time with 
the unaveraged magnitudes |ico,o(^)| = |'?^o,o(^)| in blue and 
\wi^o{t)\ = |r^i,o(t)| in green; these are propagated according 
to 0 - The overlying red line plots correspond to the mag¬ 
nitudes |ico,o(t)| = |'Co,o(t)|, = \vi^o(t)\ propagated 

according to the averaged equation (^. From the perspective 
displayed in the upper left panel the averaged and the un¬ 
averaged amplitudes seem to lie on top of one another. The 
top right panel here shows a small portion of the same plot 
zoomed in so as to show how the averaged solution in red in 
fact averages the non-autonomous equation solution (over a 
period of the nonlinearity prefactor variation) in blue. 


FIG. 8 : (color online) These plots are similar to those in Fig¬ 
ure l^but here for the double site out-of-phase compacton solu¬ 
tion; also the bottom right panel here corresponds to t = 180. 
The values of e, Q, 71 , a are the same as those seen in Figure 
1^ for the double site out-of-phase case. These values again 
give Ao ^ 1.5507 for j = 1 by 0 . The instability strength 
for this solution is max(real(z/)) ~ 0.0298. Notice that here, 
the eventual evolution leads to a single-site excitation (see the 
amplitude evolution in the top left and the eventual profile in 
the bottom right). 
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FIG. 9: (color online) We plot max(real(z/)) in order to show 
the strength of the instabilities and regions of stability. Pa¬ 
rameter values common across all four plots are L = 11 , 
e = 0.01 and Q = 71 = 10 so that a = 1. The left-most 
plot is associated to a four site symmetric solution using the 
first zero of the Bessel function with j = 1 , the middle plot 
corresponds to a four site asymmetric solution with j = 1, 
and the right-most plot is for the five site solution with j = 5. 
The axis labels are the same as those established in Figure 


case in the regular DNLS, these instability intervals end 
at a finite value of the coupling). On the other hand, for 
the out-of-phase solutions, the instability intervals arise 
for an intermediate range of the couplings, a feature to 
which we will return below. Interestingly, we also find 
that an increase in the index j of the zero of the Bessel 
function Jq decreases the size of the region of instability. 
Additional analysis has shown that for the double site 


in-phase, an increase in the parameter a (i.e. decrease 
in frequency fl/27r of the 7 function) decreases the size 
of the region of instability when plotted versus the cou¬ 
pling K, and the dc nonlinearity strength 70 as in Figure 
1^ For the double site out-of-phase, an increase in the 
parameter a (i.e. decrease in frequency 0/27r of the 7 
function) expands the region of instability when plotted 




















































versus /^, 7 o as in Figure 

The nature of the double site instabilities is more 
clearly demonstrated in Figures HI in For the double 
site in-phase solutions, unstable eigenvalues constitute a 
pair of equal in magnitude and opposite real numbers; 
see Figure where more details of the transition be¬ 
tween stable and unstable are discussed. A key feature 
of the double site in-phase (in)stability is that for fixed 
7o, small nonzero hz values give unstable solutions and 
higher k values stabilize the solution, a feature that is 
typically absent in the standard DNLS model. For the 
double site out-of-phase solutions, unstable eigenvalues 
appear as a quartet (due to the Hamiltonian nature of 
the system), i.e. two complex conjugate pairs with equal 
magnitude plus/minus real parts. This arises from the 
collision of the two imaginary eigenvalues (in this case), 
stemming from the origin with the continuous spectrum, 
a feature that is common in such out-of-phase focusing 
nonlinearity settings; see e.g. [ 22 ]. More details of the 
transitions between stable and unstable solutions for the 
out-of-phase solutions are discussed in Figure A key 
feature of the double site out-of-phase (in)stability is that 
for fixed 70 , both small nonzero and large positive k val¬ 
ues give stable solutions with the unstable region lying 
in a finite n interval that lies away from /^ = 0. For 
large k, the formerly unstable quartet grows (in imagi¬ 
nary part) faster than the band of the continuous spec¬ 
trum and eventually returns to the imaginary axis, lead¬ 
ing to spectral stability. 

For stable solutions in the averaged model the ampli¬ 
tudes of the excited sites remain steady at the predicted 
value Aq in ^ or oscillate around it, if perturbed; this 
is the same feature described above for stable single site 
and vortex type solutions. For stable solutions in the 
original non-autonomous model the amplitudes oscillate 
about the value Aq as expected. 

In Figurethe green dots show select parameter values 
for which the dynamics of unstable double site solutions 
will be explored below. Recall that in the case of the 
double-site solutions, equation © gives the amplitudes 
of the two excited sites r^o, 05 '^i,o at t = 0. In Figures]^ 

plots of the evolution of the amplitudes of the unsta¬ 
ble double site solutions over the propagation parameter 
t are shown in detail, according to the averaged equa¬ 
tion and the original non-autonomous equation §■ 
The unstable double site solutions (both in-phase and 
out-of-phase) transition to a non-compacton state with 
oscillating phase so that the initial phase profile is not 
preserved over time. For the double site in-phase solution 
in Figure the two initially excited sites remain of order 
one amplitude as t increases (although the solution mass 
is asymmetrically distributed between them due to their 
different amplitudes after the instability manifestation). 
For sites in the vicinity of the original compact support, 
there is gain in small amplitudes increasing the footprint 
of the solution. For the double site out-of-phase solution 
in Figure as t increases the amplitude at the two center 
sites first oscillates and then one of the sites’ amplitudes 


drops down towards zero; the resulting configuration is a 
non-compacton solution with order one magnitude only 
at one site, i.e., the waveform degenerates towards a fun¬ 
damental, single-site solution. 

In Figure [^ we show plots of the dependence of the 
instability strength (i.e., of the maximal growth rate of 
the potential unstable modes) on the parameter grid for 
the four and five site solutions. The stability properties 
of both the four site symmetric (in-phase) compacton 
and the five site compacton configurations follow a simi¬ 
lar pattern to that of the double site in-phase solutions; 
compare the left-most and right-most plots of Figure to 
the left column of Figure That is, for fixed 70 as k in¬ 
creases from zero the solutions are immediately unstable 
with real eigenvalues and then they eventually stabilize 
for higher values of the coupling n. The eigenvalues for 
the unstable solutions appear in real pairs similar to that 
which is seen in Fig. [^ While the double site in-phase 
solutions have only one pair of real nonzero eigenvalues, 
the four site symmetric solutions have three pairs of real 
nonzero eigenvalues and the five site unstable solutions 
have four pairs. This is in line with the expectations from 
the standard DNLS model [ 22 ] . 

The stability properties of the four site anti-symmetric 
solutions (with adjacent sites being out-of-phase) follow a 
pattern similar to that of the double site out-of-phase so¬ 
lutions, where for small enough k the solutions are stable. 
Increasing n one finds a finite (intermediate coupling) in¬ 
terval of instability, and then for higher n the solution 
stabilizes again. To see the similarity between the four 
site anti-symmetric solution stability and the double site 
out-of-phase stability compare the center plot of Figure 
l^to the right column of Figure Plots of the eigenval¬ 
ues in the complex plane look similar to the plots in Fig. 
[^ While the double site out-of-phase solutions feature 
for small coupling a single imaginary eigenvalue pair and 
thus have only one potential quartet of unstable eigen¬ 
values, the four site anti-symmetric solutions have three 
imaginary pairs for small k and may eventually feature 
up to three eigenvalue quartets, again in line with what 
one may expect in the standard DNLS case [ 22 ] . 


V. CONCLUSIONS FUTURE CHALLENGES 

Summarizing the findings of the present work, we have 
illustrated through a combination of analytical consid¬ 
erations and numerical results the existence of 2D com- 
pactons of the discrete nonlinear Schrodinger equation in 
the presence of fast periodic time modulations of the non¬ 
linearity. In particular, we showed that single site multi¬ 
dimensional compactons are very robust excitations, two- 
site stationary compactons, both of symmetric and anti¬ 
symmetric type, are also quite generic and can be sta¬ 
ble in a wide region of the parameter space (of the cou¬ 
pling and nonlinearity prefactors). Four-site compactons 
have been found to be always stable only in the vor¬ 
tex state, an unusual feature which is fundamentally dis- 
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tinct from the case of the standard discrete nonlinear 
Schrodinger model. The five site compactons may also 
feature instabilities but can be controllably again sta¬ 
bilized for suitable parametric intervals of the tunneling 
constant and/or dc-nonlinearity strength. These findings 
were not only obtained for the effective averaged DNLS 
equation that was derived herein, but were also confirmed 
by means of direct numerical simulations of the original 
non-autonmous DNLS model. 

It is relevant to point out that excitations of the 
form considered herein should in principle be accessi¬ 
ble to state-of-the-art current experiments in the optical 
(waveguide array) or atomic (Bose-Einstein condensate) 
realm. For instance, in the BEC context, such states 
could be obtained by considering a 2D array of ^^Rb or 
^Li created with a deep two-dimensional optical lattice. 
In this case, the corresponding prototypical mean-field 
model for the bosonic wavefunction 7 / would read 

ihi>t = - — V‘^'ip + g2D\'ip\'^'ip + Voism^{kx)+sm^{ky))'ip, 
Zm 

(19) 

where g 2 D is the effective nonlinear prefactor; see, 
e.g., [27]. Vb characterizes the strength and k the period 
of the optical lattice. In the (superfluid) limit of a deep 
optical lattice, e.g. Vb ^ Er = k?k^ j2m^ this equation 
reduces to the DNLS model |23]. For ^Li one could use 
the Feshbach resonance technique at the external mag¬ 
netic field of Bq = 738G with the width A5 = — 170G 
[28] to modulate the interactions. The background scat¬ 
tering length can be taken as ^ — 20aB, where as 
denotes the Bohr radius. The time dependence of the 
scattering length as{t) = Uso + Ugi cos(ci;t) follows from 
the variation of the magnetic field in time according to 

as{B)=abg (1- • 

For B{t) = (970 ± 155 cos(co’t))G, we have Uso = 
—baB^dsi = lOa^. Using suitable lattice parameters 
such as 27r//c = 0.5/im and Vb = 15Fi^, we can ensure be¬ 
ing in the regime of applicability of the DNLS, enabling. 


under the above conditions, the experimental realizabil¬ 
ity of the higher-dimensional discrete compactons. 

On the other hand, there are numerous interesting 
themes for future theoretical investigations. It would be 
interesting to examine in more detail the origin of fea¬ 
tures that are fundamentally different between the aver¬ 
aged model considered herein and the standard DNLS 
model. These include among others (a) the restabi¬ 
lization of in-phase configurations; (b) the absence of 
instabilities for the discrete vortex configuration. Ad¬ 
ditionally, it would be especially interesting to explore 
the approach of the model to the continuum limit of 
large coupling. It is well-known that collapse features 
emerge as the DNLS model approaches the continuum 
limit [22|, although this may happen in unconventional 
ways for non-standard discretizations of the model; see 
e.g. [29]. It would be intriguing to explore the prop¬ 
erties of the present discretization (and of the original 
non-autonomous model) as this limit is approached. Fi¬ 
nally, it would also be of interest and relevance to explore 
three-dimensional configurations in analogy to ones of the 
standard DNLS [22] and to examine their stability prop¬ 
erties. Such studies are currently in progress and will be 
reported in future publications. 
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